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Abstract 

A Costas latin square of order n is a set of n disjoint Costas 
arrays of the same order. Costas latin squares are studied here from 
a construction as well as a classification point of view. A complete 
classification is carried out up to order 27. In this range, we verify the 
conjecture that there is no Costas latin square for any odd order n > 
3. Various other related combinatorial structures are also considered, 
including near Costas latin squares (which are certain packings of 
near Costas arrays) and Vatican Costas squares. 
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1 Introduction 



A Costas array of order n (or side n) is an n x n array of dots and empty 
cells such that 



1. there are n dots and n(n — 1) empty cells, with exactly one dot in 
each row and column, and 

2. all the segments between pairs of dots differ in length or in slope. 

Costas arrays were introduced by J. P. Costas; see [TJ HI [7] for early results 
and historical remarks. Extensive surveys of Costas arrays can be found in 

lain]. 

Two Costas arrays of order n are disjoint if there is no cell in which 
both arrays have a dot. Let D(n) denote the maximum number of mutually 
disjoint Costas arrays of order n. The main topic of this paper is the study 
of D(n). Obviously D(n) < n 2 /n — n. We are particularly interested in 
cases where this upper bound on D{n) is attained; in those cases, we say 
that we have a full set of mutually disjoint Costas arrays. 

Example 1.1. We present all the Costas arrays of order 4. The four 
arrays in the first row are mutually disjoint. 



For notational convenience, Costas arrays are often presented using a 
certain one-line notation. Given a Costas array of order n, let n(i) = j 
whenever the array contains a dot in cell (i,j)- A Costas array of order n 
can be presented as the permutation ir = (7r(l), 7t(2), • ■ ■ , 7r(n — l),7r(n)). 



2 



We call this the permutation representation of a Costas array. Using this 
notation, the Costas arrays in the first row ofExample ll.il have permutation 
representations (1,2, 4, 3), (4, 3, 1, 2), (2, 1, 3, 4) and (3, 4, 2, 1), respectively. 

A Costas latin square of order n, denoted CLS(rt), is a latin square of 
order n such that for each symbol i, 1 < i < n, a Costas array results 
if a dot is placed in the cells containing symbol i. Clearly a CLS(n) is 
equivalent to n disjoint Costas arrays of order n. Costas latin squares were 
first defined and studied by Etzion [3J. 

Example 1.2. The first four disjoint Costas arrays of order 4 given in 
Examvle \l.l\ are mutually disjoint and lead to the following CLS(4). 
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3 



Etzion [3] defines a near Costas array of order n > 2 to be an n x n 
array of dots and empty cells such that 

1. there are n — 1 dots and n 2 — n + 1 empty cells, with at most one dot 
in each row and column, and 

2. all the segments between pairs of dots differ in length or in slope. 

As with Costas arrays, we will say that two near Costas arrays of order n 
are disjoint if there is no cell in which both arrays have a dot. 

Let D near {n) denote the maximum number of near Costas arrays of 
order n. Observe that D near (n) < L^ryJ = n + 1 for all n > 3. A set of 
n + 1 mutually disjoint near Costas arrays of order n > 3 is called a full 
set of arrays. In this case, if the filled cells of the ith near Costas array is 
given the symbol i, then superimposing all these near Costas arrays yields 
an n x n array on the symbol set S = {1, 2, . . . , n + 1} with exactly one 
empty cell. Such a square is termed a near Costas latin square of order 
n, and abbreviated as an NCLS(n). A somewhat stronger definition could 
also be given, where one of the n + 1 arrays is a Costas array and there is 
no empty cell. 
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Example 1.3. A near Costas latin square of order 6 (a NCLS(6)). 
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The roles of the rows, columns and symbols of a latin square can be 
interchanged. By interchanging the roles of rows and symbols, a Costas 
latin square of order n can be transformed into an (equivalent) latin square 
L = (hj) with the property that for any 1 < k < n — 1, hj+k ~ h,j = 
k,m+k — k,m implies that j — m. We call such a square a row-Costas latin 
square. Analogously, interchanging the roles of columns and symbols, we 
get the definition of column- Costas latin squares. 

We define one final related concept. A Vatican square of order n is 
a latin square L — (hj) of order n with the property that for any fixed 
d, 1 < d < n — 1, all ordered pairs of the form (h^,h,j+d)) 1 < * < rl i 

1 < j ' < i ~ ^ are distinct. A Vatican Costas square is a latin square 
that is both a Vatican square and a Costas latin square. Vatican squares 
were introduced in [§j while Vatican Costas squares were first considered 
by Etzion in [3]. In a Vatican Costas square one can interchange the roles 
of the rows, columns, and symbols to get other equivalent structures. 

The paper is organized as follows. In Section [2] some basic constructions 
for Costas arrays are surveyed and applied to the problem of constructing 
Costas latin squares; near Costas arrays and near Costas latin squares 
are also considered. A computational study is carried out in Section [3] to 
enumerate Costas latin squares of small order; also other related structures 
defined above are considered. 

2 Constructions 

2.1 Costas Latin Squares 

One way to construct Costas latin squares would seem to be to start with 
an initial Costas array of order n and then move each of the dots to the 
right k cells (mod n) to construct the fcth array (where < fc < n — 1). 
Certainly these arrays will be mutually disjoint, but in general they will 
not be Costas arrays. However, this technique will succeed provided that 
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we start with a singly periodic Costas array, which we now define. 

Let annxn array be left-right extended by putting a dot in cell (i, j+n) 
if and only if there is a dot in cell (i, j). A singly periodic Costas array of 
order n is an n x n array such that every n x n window in its left-right 
extension is a Costas array of order n. The following result is now obvious. 

Theorem 2.1. If there exists a singly periodic Costas array of order n, 
then there exists a CLS(n). 

Next, we describe the well known Welch construction for Costas arrays. 
Let p be prime and a be a primitive element in the field F p . Let n = p— 1. 
A Costas array of order n is obtained by placing a dot at (i, j) if and only if 
i = a? , for < j < n, and i = 1, . . . , n. We get the following result, which 
is implicit in [5] and explicit in [3J. 

Corollary 2.1. There exists a CLS(p — 1) for all primes p. 

Proof. It is known that Costas arrays constructed via the Welch construc- 
tion are singly periodic. Apply Theorem 12. II □ 

Remark. In [5], this method was used to construct a Vatican square. 

Using 2 for the primitive element in F5, one obtains the following Costas 
array of order 4, given on the left below. (Note here that the columns are 
labeled 0, 1, 2, 3, while the rows are labeled 1, 2, 3, 4 from the bottom.) This 
Costas array is singly periodic. On the right is the CLS(4) obtained from 
this array. It is interesting to note that this Costas latin square is different 
from any dihedral transformation of the CLS(4) presented in Examplc ll.il 
(e.g., consider the front and back diagonals). 
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Unfortunately, the only known singly periodic Costas arrays of order n 
come from the Welch construction when n = p — 1. It is also known that a 
singly periodic Costas array of order n does not exist if n > 1 is odd, or if 
n e {8,14,20}. 

Corollary 12.11 gives an infinite class of CLS(n). Etzion 3 gave another 
infinite class which contains no new orders but which are inequivalent to 
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the previous class. Etzion's second infinite class is based on the Golomb 
construction, which we now describe. 

Let a and (3 be primitive elements in a finite field ¥ q . Denote n = q — 2. 
For 1 < i < n, 1 < j < n, place a dot in cell of an n x n array if and 
only if ol % + /3 J ' = 1. Denote the resulting array by G{a,f3). Then G{a,f3) 
is a Costas array of order n. 

Etzion's infinite class is constructed as follows. Suppose that q = 2 r 
and suppose that q — 1 is a (Mersenne) prime. Let a £ ¥ q be a primitive 
element. Note that a 1 is primitive for 1 < i < n because g — 1 is prime. It 
is not difficult to verify that the n Costas arrays G(a, a 1 ) are disjoint, for 
1 < i < n. This yields a CLS(2 r — 2) whenever 2 r — 1 is a Mersenne prime. 

2.2 Near Costas Latin Squares 

We begin with a simple observation. 

Lemma 2.2. Suppose there exists a CLS(n). Then D near (n) > n. 

Proof. Given a CLS(n), we have n disjoint Costas arrays of order n. Re- 
move one dot from each of these n arrays, yielding n disjoint near- Costas 
arrays of order n. □ 

Corollary 2.3. If p is a prime, then D near (p — 1) > p — 1. 

We now provide another way of obtaining n — 1 near-Costas arrays 
of order n — 1. This method works whenever n is a prime power; it is 
based on the Lempel construction which is the special case of the Golomb 
construction where /3 = a. We will use the following slight variation: Let 
q be a prime power and label the rows and columns of q — 1 square arrays 
of side q — 1 by the integers 0, . . . , q — 2. The arrays will be named as 
A c , c 6 F*. For each such array A c , place a dot in cell (i,j) whenever 
a % + a 3 — c. Clearly these q — 1 arrays are all disjoint and symmetric. Each 
array contains q — 2 dots which miss exactly one row and column (namely, 
row and column x, where a x = c). 

We now show that each of these arrays is a near-Costas array. This 
can be done by modifying the proof of [H Theorem 2]. Suppose that A c 
contains dots in cells (i, j), (i + u, j + v), («', j') and (i' + u,j' + v), where 
* 7^ J 7^ j' an d w, v ^ 0. Then 

a 1 + a 3 = a l+u + a j+v = c/ + a 3 ' = a 1 '+" + a> ' +v = c. 
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We then obtain 

a %+u + a v (c- a 1 ) = a l ' +u + a v (c~ a 1 ') = c, 

and therefore 

a i+u _ a i+v = a i'+u _ a i'+v = c ^ _ a vy 

We have c ^ 0, and a v ^ 1 follows because v ^ 0. Hence, 

a\a u ~a v ) = a 1 ' {a 11 - a v ) ^ 0. 

Notice that u ^ v because a 1 (a" — a 11 ) 7^ 0. However, since u ^ v, we 
obtain a 1 = a 1 , which is a contradiction because i ^= i . 

Summarising the above discussion, we have the following theorem. 
Theorem 2.4. If q is a prime power, then D near (q — 1) > q — 1. 

Remark. Theorem 12.41 is more general than Corollarv l2.3l in that it applies 
to prime powers rather than just primes. 

Example 2.1. Let q = 7 and take a — 3. We present the superposition of 
Ai,...,A 6 : 
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Next, we use the Golomb construction to construct disjoint near Costas 
arrays of order q — 2, for certain prime powers q = 3 mod 4. Since q = 
3 mod 4 is a prime power, —1 is a quadratic non-residue modulo q. Suppose 
also that q = 2p + 1 where p is an odd prime. Then ¥ q has (j){2p) = p — 1 
primitive roots, say a^, fe = 1, . . . ,p— 1. 

Denote 7 = a^, <5 = and e = a m , where £ ^ m. We will show 
that G(7, (5) and G(j, e) contain exactly one common dot. Suppose to the 
contrary that 7* + = j l + e J = 1. Then $ = e 3 '. Denote £ = £/e; then 

= 1. It is clear that ord(£) | 2p, so ord(£) e {l,2,p, 2p}. We consider 
each possibility in turn. 

case 1: ord(£) = 1. Here £ = 1, so S = e, which is a contradiction. 
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case 2: ord(£) = 2. Here £ = — 1 and S = — e. However, —1 is a quadratic 
non-residue, so S and e cannot both be primitive elements. This is a 
contradiction. 

case 3: ord(£) = p. Here C is a quadratic residue and # = e> = —1. This 
implies that j — p and i is the unique clement such that 7' = 2. 

case 4: ord(£) = 2p. Here ( is a primitive element and <5 J = e 3 = 1. This 
implies that j = 2p, which is a contradiction because j < 2p — 1. 

It follows from this discussion that the p— 1 Costas arrays G(ai,ai) 
(1 < i < j? — 1) all contain a single common dot. On removing this dot, we 
obtain p — 1 disjoint near-Costas arrays of order 2p — 1. Summarising, we 
have the following. 

Theorem 2.5. Suppose that q = 3 mod A is a prime power and suppose 
that (q — l)/2 is prime. Then D near (q — 2) > (q — 3)/2. 

Theorem 12.51 is illustrated by the following example. 

Example 2.2. Let q = 11, so n = 9 and p — 5. The four primitive 
elements modulo 11 are 2, 8, 7 and 6. W^e present the superposition of 
G{2,2), G(2,8), G(2,7), G(2,6). TTie unique common dot is in position 
(1,5). If this dot is deleted, then we get four disjoint near Costas arrays of 
order 9. 
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Remark. It may be the case that various Costas arrays of the form G(a, a) 
are disjoint. When q = 11, one finds that the four Costas arrays G(2,2), 
G(8,8), G(7, 7) and G(6,6) are disjoint. The resulting array is presented 
in the next example. 
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Example 2.3. The four disjoint Costas arrays G(2,2), G(8, 8), G(7, 7) 
and G(6, 6). 
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It is not obvious how to derive a simple algebraic formula to determine 
when Costas arrays of the form G(a, a) are disjoint. So at present we have 
no theorem that generalises the previous example. 

3 Classification of Squares 

The concept of equivalence (or isomorphism) is central in any classification 
of combinatorial structures. Two Costas arrays are said to be equivalent if 
there is a mapping in the dihedral group of rotations and reflections of the 
square (which has order 8) that maps one array onto the other. Two Costas 
latin squares are equivalent if there is a mapping in the same group that 
maps every Costas array of one square onto a Costas array of the other. 
(Note that with this definition, the numbers given to the Costas arrays in 
the Costas latin square are irrelevant.) 

The set of all mappings that map a Costas array (or Costas latin square) 
onto itself forms the automorphism group of the respective structure. The 
definitions of equivalence and automorphism groups for near Costas arrays 
and near Costas latin squares are analogous. 

3.1 Costas Latin Squares 

Let us first have a brief look at Costas latin squares with very small pa- 
rameters. It is trivial to construct CLS(l) and CLS(2); these are unique 
up to equivalence. There is no CLS(3) since no Costas array of order 3 
can contain a dot in the middle cell. However, it is easy to construct two 
disjoint Costas arrays of order 3, so it follows that D(3) = 2. 

A CLS(4) is presented in Example II .21 In the next case, n = 5, it turns 
out that there is no Costas latin square of order 5. This can still be done by 
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hand (for example, with a case- by-case argument). The following example 
shows that there are four disjoint Costas arrays of order 5, so it follows that 
D(5) = 4. 

Example 3.1. Here are four disjoint Costas arrays of order 5, given in 
permutation representation: 

(1,4,5,3,2), (2,5,3,4,1), (3,1,2,5,4) and (4,3,1,2,5). 

It is interesting to note that a 5th Costas array which is disjoint to the other 
four arrays in three cells also exists. This Costas array has permutation 
representation 

(5,2,4,3,1). 

If we superimpose these five Costas arrays, we get the following array. 
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We know of no better way than a computer search to handle orders 
greater than 5. For a given order n, the largest possible number of mutu- 
ally disjoint Costas arrays can be found by solving the following instance 
of the maximum clique problem. Form a graph G with one vertex for each 
Costas array and with an edge between two vertices exactly when the cor- 
responding arrays are disjoint. Any clique in G corresponds to a set of 
mutually disjoint Costas arrays, so we may now search for the size of a 
maximum clique. The Cliquer software [10) was used to search for (max- 
imum) cliques in the current work. All Costas arrays of order between 1 
and 27 can be obtained in electronic form from |llj . 

Whenever there are n disjoint Costas arrays of order n, one can find all 
cliques of size n in the above mentioned graph to get all possible Costas latin 
squares. However, utilizing the fact that the arrays then cover all the cells 
of the square, such cases are for performance reasons preferably considered 
as instances of the exact cover problem. In the exact cover problem, we are 
given a set S and a collection U of its subsets, and the task is to form a 
partition of S by using sets in U. Here, S is the set of all cells of the square, 
and U has one set for each Costas array (consisting of the cells where there 
are dots). The libexact library [9] was used here to solve instances of the 
exact cover problem. 
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The computational results for all orders up to 27 are summarized in 
Table [TJ The columns in the table are the order n, the total number of 
Costas arrays Nc, the number of equivalence classes of Costas arrays Nee, 
the largest number of disjoint Costas arrays D(n), and, if a full set of 
disjoint Costas arrays exists, the total number of Costas latin squares Nl 
and the number of equivalence classes of Costas latin squares Nle. The 
entries for Nc and Nee have been taken from (2J [6l [12] . 



Table 1: Computational results for Costas arrays 
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Nee 


D(n) 
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Nle 
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7852 
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12 
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3761 


13 


12828 


1616 


< 12 






14 


17252 


2168 


< 13 






15 


19612 


2467 


< 14 






16 


21104 


2648 


16 


32768 


8256 


17 


18276 


2294 


< 16 






18 


15096 


1892 


18 


5832 


756 


19 


10240 


1283 


< 18 






20 


6464 


810 


< 19 






21 


3536 


446 


11 






22 


2052 


259 


22 


200 


30 


23 


872 


114 


9 






24 


200 


25 


8 






25 


88 


12 


5 






26 


56 
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27 


204 


29 


8 







It is easy to determine the automorphism groups of the classified Costas 
latin squares. The results were then validated by applying the orbit- 
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stabilizer theorem to get the total number, which coincides with the number 
in the column Nl of Table [T] 

Note that there exists a CLS(8); this was also observed in [3J. This is 
the only order n for which there exists a CLS(n) but n + 1 is not a prime 
number. An example of a CLS(8) can be found in the Appendix. 

Table Q] lends additional numerical evidence to the truth of the con- 
jecture of Etzion [3j Conjecture 2] that there is no CLS(n) for any odd 
n>3. 

3.2 Disjoint Near Costas Arrays 

For near Costas arrays, one can carry out a similar computational study as 
for Costas arrays in the previous subsection. In this case exact cover cannot 
be used in a direct way to find full sets of disjoint arrays. One possibility 
is to use a clique search to find full sets, another is to use the framework 
of exact cover for all n 2 possibilities for the empty cell in the near Costas 
latin square. 

In any case, the number of near Costas latin squares grows very quickly; 
already for order 5 this number is 978982. Consequently, we will not be able 
to present an extensive table like that in TableQ] Instead, in the Appendix, 
we give just an example of a near Costas latin square for each order up to 
n = 8. 

3.3 Vatican Costas Squares 

For some of the structures defined in the introduction, the mappings in 
the definitions of equivalence must map rows to rows (resp. columns to 
columns), that is, the mapping is in a dihedral group of order 4 rather than 
8. Vatican Costas squares, which we consider here, are such structures (so 
this specification should be added to the definition on [3j p. 149]). 

We consider one representative from each equivalence class of Costas 
latin squares classified in the current work (cf. column Nle in Table [1]), and 
check whether it is a Vatican Costas square. This check is also carried out 
for the square obtained by a 90 degree rotation. For the orders up to 27 for 
which Costas latin squares exist, that is, n = 2, 4, 6, 8, 10, 12, 16, 18 and 
22, the number of equivalence classes of Vatican Costas squares is 1, 1, 3, 
0, 2, 2, 4, 3 and 5, respectively. One square from each equivalence class is 
given in the Appendix. Etzion [3] already classified such squares of order 6. 
It is interesting to note that many of these squares have the property that 
all of the columns are translates of a single column. Also note that n = 8 
is the only order for which there is a Costas latin square but no Vatican 
Costas square (coincidentally, perhaps, it is also the only order n where 
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there is a Costas latin square but n + 1 is not a prime) . 
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Appendix 

A.l A Costas latin square of side 8 
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A. 2 Near Costas latin squares of orders 2 through 8 
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